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Quantum to classical transition and entanglement sudden death in Gaussian states 

under local heat bath dynamics 

Sandeep K Goyajj and Sibasish GhosllJ 

Center for Quantum Sciences, The Institute of Mathematical Sciences, 

CIT campus, Taramani, Chennai 600 113, India 

Entanglement sudden death in spatially separated two-mode Gaussian states coupled to local 
thermal and squeezed thermal baths is studied by mapping the problem to that of the quantum- 
to-classical transition. Using Simon's criterion concerning the characterisation of classicality in 
Gaussian states, the time to ESD is calculated by analysing the covariance matrices of the system. 
The results for the two-mode system at T = and T > for the two types of bath states are 
generalised to n-modes, and are shown to be similar in nature to the results for the general discrete 
n-qubit system. 



I. INTRODUCTION 

Entanglement is one of the basic features that distin- 
guish quantum systems from their classical counterparts 
and is the most useful resource in Quantum Information 
Theory (QIT)[l|. It is indispensible for essential quan- 
tum information tasks like quantum computation and 
teleportation as well as super-dense coding and one-way 
communication, to name a few. However, more often 
than not, quantum systems do not operate while com- 
pletely insulated from its environment. Unfortunately, 
quantum systems are inherently fragile and extremely 
sensitive to environmental interactions. It is a common 
belief that the evolution of a quantum system in contact 
with a dissipative environment results in an asymptotic 
transition to classicality |2-[J1 and thus an inevitable loss 
of entanglement, since classicality subsumes disentangle- 
ment. Recently however it has been shown p-Q that 
there exists a certain class of two-qubit states which dis- 
play a finite entanglement decay time. This phenomenon 
is aptly called Entanglement Sudden Death (ESD) and 
cannot be predicted from quantum decoherence which is 
an asymptotic phenomenon. 

ESD has important and obvious implications for the 
success of quantum information tasks. Much research has 
been done on ESD in discrete quantum states. Qasimi 
et al. 3 have shown that ESD occurs for a class of two- 
qubit states called 'X' states at non-zero temperatures. 
The authors of the present paper have shown that there 
is ESD for all n-qubit states at finite temperature [9| and 
have also given a sufficiency condition for the existence of 
ESD for a system consisting of n number of d-dimensional 
subsystems as well. The study of the evolution of entan- 
glement as well as that of its sudden death in continuous 
variable systems has also received a great attention, for 
example, the paper by Dodd et al. [l^l deals with ESD 
from the point of view of a separable representation of 
the joint Wigner function of two-mode Gaussian states. 
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Their results do not have an explicit presence of tempera- 
ture. Also, the time to ESD for the maximally entangled 
EPR state is shown to be a lot smaller than that of any 
other initially entangled two-mode Gaussian state. This 
is counter-intuitive as one would expect the ESD time for 
a maximally entangled state to be the largest. In [ll| . 
Diosi has given a bound on the time of ESD by using 
a theorem on 'entanglement breaking quantum channels 
p^ - [l4l |'. In [i^ Marek et al. have addressed a different 
problem. In this paper they obtained a class of states 
which is tolerant against the decoherence at zero tem- 
perature. People have studied the decoherence in infinite 
dimensional systems interacting with different kinds of 
bath and with different system-environment models |16l - 

Since classicality subsumes disentanglement, the ESD 
problem can be embedded into the corresponding 
quantum-to-classical transition problem. Thus, the time 
taken for the system to attain classicality will be an up- 
per bound on the time at which entanglement dies. We 
use the criterion due to Simon [20|, l21| for checking the 
classicality of Gaussian states in terms of their covari- 
ance matrices V , which enables us to characterise the 
quantum-to-classical transition. If I be the covariance 
matrix corresponding to the vacuum state, then Simon's 
criterion states that Gaussian states attain classicality if 
and only if the difference V — 1 becomes positive semi- 
definite. Which is equivalent to Sudarshan-Glauber P 
function being positive. We use this condition to find 
out the transition time to classicality (tc) and hence the 
time to ESD. Diosi et al. [221 discusses about the positiv- 
ity of Wigner function and P function under Markovian 
evolution. We find that the behaviour of the Gaussian 
system, as regards to ESD, depends on the state and the 
temperature of the bath it is coupled to. We have shown 
that all Gaussian states show ESD at all temperatures 
T > 0, whereas some do not at T = when coupled to 
a thermal bath. However, all states again show ESD at 
T = when the state of the bath is squeezed thermal. 

We begin this article by writing down the master equa- 
tion for a single harmonic oscillator (Sec. |ll| in con- 
tact with a thermal bath of infinitely many oscillators, 
taken at some finite non-zero temperature. We express 



the states of this system in the Sudarshan- Glauber P- 
representation and derive the corresponding quantum 
analog of the classical linear Fokker-Planck equation 
(Sec. imp . This enables us to use covariance matrices 
to characterize classicality of the single-mode system in 
the presence of thermal baths (Sec. IIV|) . Once we have 
this mechanism in place, wc use it to determine the action 
of the thermal bath (Sec. |V]) and squeezed thermal bath 
(Sec. IVl)) on a two-mode Gaussian state. Using Simon's 
criterion, we calculate the time to classicality and hence 
are able to determine when a system will show BSD. Fur- 
thermore, using the two-mode analysis, we are able to 
predict the possibility of ESD in any n-mode Gaussian 
state since we make the important observation that Si- 
mon's criterion y — I > is independent of the number 
of modes of the system. Hence we are able to generalize 
our earlier statements to n-modc systems as well (Sec. 
IVII| . Section FVlIII provides the conclusion of the paper. 



II. MASTER EQUATION FOR SINGLE MODE 

SYSTEM 



The irreversible time evolution of the state ps {t) of a 
single harmonic oscillator coupled to a bath is described 
by the Lindblad master equation (ME) |23l |. 



d 



^/3s(0 = -iuJoia' a, psit)] 

+ joiN + 1) I ap,it)a'' - -a'^apsit) - -ps{t)a'^a I 

+ 7o(^) \ a'Ps{t)a - ^aa^Psit) - -ps{t)aa) \ . 

(1) 

The initial state of the entire system is of the form 
PsiO) ® pth, where ps{0) is the initial state of the sys- 
tem and the thermal state pth is the initial state of bath. 
Since the bath we are considering is very large, the state 



of the bath will not change during the evolution. The 
first term of the ME describes free evolution generated 
by the system Hamiltonian H = ujQa^a while the rest 
are interaction terms with the bath, where the dissipa- 
tive coupling is provided through the damping rate 70. 
Here N = (e*^"" — 1)^^ is the mean number of quanta 
in a mode with frequency luq. This ME can be used to 
describe, for example, the damping of an electromagnetic 
field inside a cavity where a and a^ denote the creation 
and annihilation operators of the cavity mode. The mode 
outside the cavity plays the role of the environment with 
dissipative coupling rate 70. 



III. 



COHERENT STATE REPRESENTATION 



We can transform the master equation described by ([T]) 
into a continuous diffusion process by using the coherent 
state representation. 



Ps{t)^ / (faP{a,a*,t)\a){a\, 



(2) 

where P{a,a* ,t) is the Sudarshan-Glauber P-function 
J24L |25| and the integration in ([2]) is over the entire com- 
plex plain. The quasiprobability distribution, called thus 
because P{a,a* ,t) can take negative values for some a, 
satisfies the normalization 

tTsPs{t)= f d^aP{a,a*,t) = l. (3) 

Substituting ^ into ([T|) and using the properties 



a|a)(a| = a|a)(Q;|, 



(4) 



we get the following equation for the evolution of 
P{a,a*,t): 



J 



d_ 
di 



P{a, a* ,t) 
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da* 



P{a,a*,t)+-foN 
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dada' 



-P{a,a\t). 



(5) 



This is structurally similar to the classical linear Fokker- 
Planck equation [26| and can be solved using the Gaus- 
sian ansatz 



P{a,a\t) 



1 



7r(T2(f) 



exp 



a- 



m 



21 



a^{t) 



(6) 



given the initial condition P{a,a* ,0) = 5^{a ~ ao), 
where |ao) is the initial coherent state. Here the 
mean amplitude is given by /3(i) = j d'^aaP{a,a* ,t) = 



Q,pg(-i"o-7o/2)t g^-^^i ^Yi^ variance is (j'^{t) = 7V(1 -e^^"*). 
Since we arc interested in dissipation, we choose to ig- 
nore the bjQ term which only contributes to free evolu- 
tion. The linearity of the evolution map ensures that any 
density matrix p(0) = / P(A, A*,0)|A)(A|d2A will evolve 
to p{t) = j P{a,a* ,t)\a){a\d^a, where P(a, a*,i) = 

/ P(A, A* , 0) cxp [~ '°;fg'' ) d^P with p = Ae-^°*/2. 
Therefore, thermal evolution of the P-distribution mani- 



fests as a convolution of the P and the thermal distribu- 
tion Pth = exp(-|a|Va2(^))^ 

IV. COVARIANCE MATRIX 

Our aim is to express the evolution of the state in terms 
of the covariance matrix, and to that end we write down 
the symmetric characteristic function x(a, t) [23, 12^ 

X{a,t) ^ e-^ T[P {a, a* ,t)] ^ cxp ( j , (7) 

where J" denotes the Fourier transform and V is the 
covariance matrix. The vector X^ = [q, p) where 
q = -j=(^a + a*) and p = -7^{oi* — en) are the position 

and momentum variables. The evolution of x(q^) t) fol" 
lows from that of P{a, a* , t) (equation[5]). We have, from 



X(a,i)e— -x(/3e-^*,0)e-^^e-^r^. (8) 



where A 



-t«*/2L and B 



The covariance matrix V can be written as: 



V,, 



{XX^),jT-\{a,t)dPqd'p^ 



(9) 



and hence 



V{t) = e-^°V(0) + (y + l) (1 - e"^°')ll- (10) 

A Gaussian channel is a map that takes Gaussian 
states to Gaussian states, an example of which is the evo- 
lution, given by Eq. ([T]) for a harmonic oscillator coupled 
with a thermal bath. The evolution of the covariance 
matrix V of the system, under the action of a general 
Gaussian channel, can be characterized by two matrices 
A and B: 



Vf = AVA^ + B, 



(11) 



where B is a positive operator [23|- For a thermal bath, 
we have from ^, A = e''^"*/^! and B = {f + l) (1 - 
g-7ot^|^ We can hence characterize the action of a ther- 
mal bath completely using these two matrices. The char- 
acterization, given by Eq. (fTT]) . guarantees that V{t) of 
Eq. ([TU]) is a bona fide covariance matrix for all finite 
time t. 



V. ESD OF TWO-MODE GAUSSIAN STATE 

Consider a two-mode system coupled to two identical 
local thermal baths. Let us assume that the initial two- 
mode (4 X 4) covariance matrix Vq represents an entangled 
state. Its subsequent evolution is given by 



V{t) = {A®A)V{0){A®A) 

'N 

y + 1 



-70* 



V{0) 



+ {B®B) 
(1 - e-''°')I. 



(12) 



'N_ 



e <"-i-\l4^ ana B = (^3- + ij U - e 
From a quantum optics point of view, we know [2C 
in view of the P-representation, that the state p = 
/ (i^a(i^/3P(a,/3)|a)(a| ® \P){P\ is classical only when 
P{a, P) is non-negative for all a and /3. This interpre- 
taion of classicality can be translated into the language 
of the covariance matrix V{t). Thus, a two-mode Gaus- 
sian state will be classical at some time t if and only if 
V{t) > I [20]. Clearly, being entangled, the initial co- 
variance matrix satisfies: 1^(0) < I. However, since the 
evolution (fT2|) is dissipative, the system will attain classi- 
cality after a time tc so that V{tc) > I. This condition is 
equivalent to nmin{tc) > 1, where n„iin{tc) is the small- 
est eigenvalue of V{tc)- Using (fT2]l . it can be written as 
the evolution of n„iin (0) which is the smallest eigenvalue 
of Vq- Since ^(0) is not classical, n,„i„(0) < 1. 



.(tc) 



.(0) 



N 



+ 1(1 



-lota 



(13) 



Here we find that nminitc) > 1 always for tc > 
— ::^ In I jv-f2-2rt — ~7oT ) ■ ^^ appropriate choice of Vq al- 
lows us to make nmm(O) arbitrarily small (which is the 
case in EPR states, which are maximally entangled ini- 
tially) and thus get an upper bound on the transition 

time tc, given by t,nax = -^^'<^ (l^j • For T > 0, we 
have A^ > and hence tmax is finite and non-zero. This 
proves that there is always ESD at non-zero tempera- 
tures. However, when T — 0, i.e. N — 0, we have 
tmax — >■ 00 and hence no quantum-to-classical transi- 
tion is seen at finite times. This does not, however, rule 
out ESD since non-classicality does not necessarily imply 
that there is entanglement. 

For T = 0, let us consider a particular form of the 
initial covariance matrix representing a symmetric two- 
mode Gaussian state. 



V{0) 



/ n k^ 

n —ky 

k^ n 

\ ~ky n 



(14) 



Here kx and ky are positive and satisfy n^ — 
{max{krc,ky})'^ > 1 (in order that V{0) is a bona fide 
covariance matrix). This covariance matrix represents 
an entangled state when {n — kx){n — ky) < 1 [30[. 
Under evolution (1121) and taking k^ — ky, we have 



i{t) 



-T»* + {1- e-T"*) and k^t) 



-'^«*. At 



the quantum-to-classical transition time tc, the states be- 
come separable: {n{tc) — kx{tc))'^ > 1- This condition is, 
in this case, equivalent to n{tc) — kx{tc) > 1- However, 
we have verified that there is no positive finite value of 
tc that satisfies this condition. Hence, we see that for 
the zero-temperature case, there are states given by p^ 
that do not show ESD. This can be easily generalized to 
the case when k^ 7^ ky. This result is similar to the two- 
qubit case where we sec that ESD occurs for all states 



at non-zero temperature. For T — 0, there exist states 
which do not show ESD Q. 



VI. SQUEEZED THERMAL BATH 



does not become zero with temperature T and hence a 
quantum-to-classical transition always happens at zero 
temperature for the squeezed thermal bath which ensures 
ESD. 



If the initial state of the bath is a squeezed thermal 
state, then the master equation ^ is replaced by |31j : 

— Ps(i) = -iujQ[a'a,ps{t)] 

+ joiN + 1) I ap,{t)a'f - -a'^apsit) - -Ps{t)o)a \ 

-t- 7oA^ < a^ ps{t)a - -aa^ ps{t) - -ps{t)aa^ 

M* Vlapa — a" p — pa \ 

-|Af{2aVt-(at)%-p(at)^}, (15) 

where M = -i sinh(2r)e*'^(27Vt,, + 1) and 2N + I = 
cosh(2r)(2A^t/i + l). The quantity Nth is the average num- 
ber of photon in the thermal state and r and (j) are the 
squeezing parameters. Repeating the same procedure as 
earlier, we finally write down the evolution in terms of 
covariance matrices: 



V{t) = e-^°^V({)) + (1 - e-~^°')Vo, 
where Voo is given by 



(16) 



Kc 



f + 1 + Re{M} Im{M} 



Im{M} 



f + 1 - Re{M} ' ■ (^^^ 



In this ease, for a squeezed thermal bath, A ~ e '^"^1 and 
5 = (1 — e~'^°*)Vrxi- If there is no squeezing {i.e. if we 
set M = and TV = Nth), we get V^ = {N/2 + 1)1 and 
thus recover the unsqueezed result namely, Eq. p2|) . 

Thus the evolution of the two- mode Gaussian state can 
be written as 



V{t) = e-T°V(0) + (1 - e-^°'){V^ © Foo). 



(18) 



If the smallest eigenvalue of V{0) is n(0) then the classi- 
eality condition V{t) > I can be written as [33 |: 

n{t) > e-'">*n{0) + (1 - e"^"*) (y + 1 - l^^l) > 1, 

(19) 

i.e,i > In -^ , 

- 70 V^ + 2-2|M|-2n(0)y 

(20) 

where in Eq. p^ . n{t) is the smallest eigenvalue of 
V{t) and ( Y + 1 ~ l-^l) i^ ^^^ smallest eigenvalue of Voo. 
Since n(0) < 1, ic is finite and positive. Therefore, the 
transition to classicality and thus ESD is ensured for the 
squeezed thermal bath for non-zero temperatures. How- 
ever, unlike the thermal bath case, we see here that N 



VII. ESD FOR n-MODE GAUSSIAN STATES 

In this section we generalize the result mentioned pre- 
viously for n-mode Gaussian states in the presence of 
local thermal and squeezed thermal baths. Consider a co- 
variance matrix Vn (0) corresponding to an n-mode Gaus- 
sian state. Let the matrices A and B characterise the 
single mode Gaussian channel for a given bath. The evo- 
lution of the covariance matrix can be written as [as a 
generalization of Eq. P^ )] , 



Vn{t) 



If n(0) is the smallest eigenvalue of Kt(0), the classical- 
ity condition Vn{t) > I gives rise to an equation which is 
same as Eq. p^ ( Eq. P^ ) representing nmin{t){n{t)) 
as the smallest eigenvalue of Vn{t) in the case of the 
thermal (squeezed thermal) bath. Therefore, all n-mode 
Gaussian states show quantum-to-classical transition at 
finite temperature. 



(A © A © • • • © A)V„{0){A © A © • ■ 


■■®Af 


+ {B®B®---®B), 


(21) 



VIII. 



CONCLUSION 



In this article we have analysed the phenomenon of en- 
tanglement sudden death for the case of Gaussian states 
coupled to local identical thermal and squeezed thermal 
baths. We have mapped this problem of non-asymptotic 
decay of entanglement to the quantum-to-classical tran- 
sition phenomenon by noting that classicality subsumes 
separability. Using the powerful criterion given by Si- 
mon, we have characterized the classicality of the Gaus- 
sian states and calculated the time taken to attain clas- 
sicality. We have shown whether ESD is possible at zero 
temperature as well as at non-zero temperatures for the 
two different types of baths considered. We have been 
able to generalize our two-mode result to a general n- 
mode system by exploiting the mode-independence of Si- 
mon's criterion. Finally, we have noted the similarity 
that the ESD results for Gaussian systems have with dis- 
crete system results Q- We believe that the same might 
hold for non- Gaussian states as well, work on which is 
currently in progress. 
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